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Abstract 

We derive dual parameterization for the DVCS amplitude to the twist-3 accuracy. 
The relation of the dual parameterization to the Abel transform tomography for 
amplitudes of hard exclusive processes is elaborated. We prove dispersion relations 
for the twist-3 DVCS amplitude in the Wandzura - Wilczek approximation and show 
that the subtraction constant is given by the D-term form factor. The same as in 
the case of twist- 2 DVCS amplitude. 
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Dual parameterization facilitates Abel tomography for DVCS amplitudes 

The dual parameterization [I] of generalized parton distributions (GPDs) [2J (for recent 
reviews of GPDs see [3, 4, 5] 6j) provides us handy and flexible tool to describe amplitudes 
of hard exclusive processes in terms of single function - so-called quintessence function. 
This function quantifies the maximum amount of information about GPDs, which can be 
extracted from data on hard exclusive processes. First phenomenological applications of 
this parameterization to DVCS give very promising results [3 [8]. In present paper we 
extend the dual parameterization to the twist-3 GPDs in the Wandzura- Wilczek approx- 
imation. 

Let us first remind basics of the dual parameterization and dwell on its relation to 
Abel tomography [9]. The leading order twist-2 amplitude of hard exclusive reactions is 
expressed in terms of the following elementary amplitudd3: 
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We see that the amplitude is given by the convolution integral in which dependence of 
GPDs on variable x is "integrated out". One can not completely restore the GPD H(x, £, t) 
from Eq. (Op) , because from this equation one can obtain the same amplitude with different 
GPDs H(x,£,t). In dual parameterization the expression for the amplitude has the form 



"As the first step we restrict ourselves to DVCS on spinless hadron. 
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The amplitude is completely determined by, so-called, GPD-quintessence function N(x, t) 
[8J and by the D-form factor D(t). The latter is : 



D{t) = Y J d n {t) = - / dz- 

n=l 



(4) 



where D(z, t) is the D-term [12] . One can check [10] that the amplitude given by Eqs. 
automatically satisfies the dispersion relation with the subtraction constant given by the 
D- form factor, as it should be on general grounds [EUdl]. The GPD-quintessence function 
N(x,t) is expressed in terms of set of forward-like function^ as: 
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N(x,t) = ^x 2v Q 2 u{x,t). (6) 

The knowledge of the LO amplitude is in one-to-one equivalence to the knowledge of 
the function N(x,t) and the D-form factor D(t), because the Eq. (j2J) can be inverted [TO] , 
i.e. the function N(x, t) can be expressed unambiguously in terms of amplitude. As the 
function N(x, t) is unambiguously restored from the amplitude it means that it contains 
maximal information about GPDs which one can, in principle, obtain from the amplitudes. 

t We call the functions Q2v{x,t) forward-like because [10J: 

• At the LO scale dependence of functions Qiv{x, t) is given by the standard DGLAP evolution 
equation, so that these functions behave as usual parton distributions under QCD evolution. 

• The function Qo(x,t) is related to the forward i-dependent quark densities q{x,t) (at t — they 
are reduced to usual parton densities q(x) measured in DIS) as: 



Qo(x,t) = q(x,t) +q(x,t) - - J (q(z,t) + q(z,t)) . (5) 

• The expansion of the GPD H(x, ^, t) around the point ^ = with fixed x to the order £ 2iy involves 
only finite number of functions Q2/j,(x,t) with [i<v. 

See original papers [H [101 EJ f° r detailed discussion of properties of these functions and their physics 
interpretation. 
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This inversion, actually, corresponds to the Abel transform tomography [9j. Let us show 
relation of Eq. (J2j) to tomography^. It is useful to make the following substitution for the 
integration variable x in Eq. (j2J): 




This substitution correspond to famous Joukowskii conformal map [H], which historically 
was used to understand some principles of aerofoil design. After this change of variables 
the expression for the imaginary part of the amplitude gets the form: 



ImA to2 (e,t) = t— M(w,t)^L=, (7) 
where the function M(w, t) is related to the function N(x, t) as: 



Af(tw,t) = JV ,t) ■ 8 

V w / y/2{l-w 2 Wl- VT^^ 

The Eq. (jZj) is typical for the Abel tomography. In two dimensions, the Abel transform 
a(y) can be interpreted as the projection of a axially symmetric function m(p) along a set 
of parallel lines of sight which are at distance y from the origin. Referring to Fig. [TJ the 
observer see the image: 

/oo 
dx m(p) . (9) 
-oo 

Now, using axial symmetry of the function m(p), we make substitution for the integration 
variable x 2 = p 2 —y 2 , as the result we obtain the following expression for the image function 
a(y): 

a(y)= Tdp 2 ^=. (10) 
Jy 2 VP - V 

This equation has exactly the form (with obvious identification of variables and functions) 
of Eq. ([3). We see that we reduce our problem about GPDs to the restoring of images 
of axially symmetric object^] from their photographs on plane. This problem is solved 
unambiguously with help of Abel transform [9]. Let us demonstrate us this directly for 
Eq. (EJ). Introducing the following notations m(w,t) = and a(£, £) = -^ImA(£,t) 

(p 2 = w and y 2 = £), we arrive at the integral equation of Abel: 

1 j... m ( w ^) 



a(£, t) = J dw 



* See Ref. [15] for the pioneering applications of the Radon tomography methods to GPDs. 
§It would be interesting to understand the deep meaning of this axial symmetry in context of GPD 
problem. Amplitudes of hard exclusive processes are photographs of some object which is axial symmetric 

in plane of variables X = J — £ and Y = 
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Figure 1: Illustration for Abel transform tomography 



which has very simple and unambiguous solution. Indeed, to solve the Abel equation we 
multiply both sides of the above equation equality by 1/V£ — z an d integrate it over £ 
from z to 1. Then changing sequence of integration over £ and w, we obtain 

Now differentiating that expression over 2, one can write 

m(w, t) = — / a£ 



g?£ ^' L — 7r I dw m(w,t). 



71 dlU J w y/£ — W 

This solves the Abel tomography problem. Making necessary trivial change of variables 
we obtain the expression for the GPD- quintessence function N(x, t) in terms of the 
amplitude [TO] : 



The beauty of the Abel transform lies in its exactness and conciseness. Many effective 
numerical methods to implement the Abel transform were developed, see e.g [PT[ . For 
instance, Chan and Lu [IB] suggested the form of Abel transform which avoids derivatives 
(which can be noisy when applied to data) of the "image function" (ImA in our case). 
Our Eq. (fTTl) in representation of Chan and Lu has the form: 
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Dual parameterization for twist-3 DVCS amplitude 

The DVCS twist-3 amplitude was computed in Refs. [HI [20] EU [22], it includes the 
contributions which are suppressed by power of 1/Q relative to the twist-2 amplitude. 
The inclusion of such terms is mandatory to ensure the electromagnetic gauge invariance 
of the DVCS amplitude to the order \/—i/Q, also these contributions can give sizable 
contributions to the DVCS observables [23]. At the order 1/Q the DVCS amplitude 
depends on a set of new GPDs. In Refs. [2H [25], [22], [27] it was shown in that in the so- 
called Wandzura-Wilczek (WW) [24] approximation these new functions can be expressed 
in terms of twist-2 GPDs. The WW approximation corresponds to neglecting the nucleon 
matrix elements of the 'genuine twist-3' operators of the type tpGip, i-e. neglecting the 
nonperturbative quark-gluon correlations in the nucleon. In principle this approximation 
is not justified by any small parameter of the theory. However, recent measurements 
of the polarized structure function gr(x) [29j show that the WW approximation works 
pretty well in this case. Also the estimates of the 'genuine twist-3' contributions polarized 
structure functions gr(x) and hi(x) [30] [31], and to twist-3 GPDs [32] in the theory of 
instanton vacuum (which is a model of non-perturbative quark-gluon correlations) showed 
that these contributions are parametrically suppressed relative to the 'kinematical' part 
of the twist-3 contributions by the packing fraction of instantons in the vacuum. 
To the twist-3 accuracy the DVCS amplitude has the form: 



T" u = — — / dx 



2P-QJ U-rr-iO £ + x - iO 
+ g^(P ■ Q) - \p^ v ± + \p u K 

[H 3 (x,0 + -Ha(x,0]A1\^P» + Q' 
x \ 



H(x, -2^P V - P»Q V - P^Q" 



(13) 



where P = (p + p')/2, Q = (q + q')/2, A = q — q' = p' — p = —2^P + Aj_, while p,p', q, q' 
are the initial and final momenta of pion and photon, respectively. The three terms in the 
square brackets are individually gauge invariant up to order Aj_, i.e. q' u T^ v = 0(A\) = 
q fJi T fJ,l/ . As P and Q are longitudinal and transverse, the first term corresponds to 
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transverse polarization of the virtual photon, whereas the second term to the longitudi- 
nally polarized virtual photon The second term contains the new twist-3 contributions to 
the DVCS amplitude, which are defined as: 



2?T 



A 



A 



(p'\^[ -^nAYlBipf )\p) 



P»H{x,Z) + A1_H a (x,Z), (14) 
^ATV^f^). (15) 



Here the light-cone vector n is normalized as n ■ P = 1 

We introduce elementary twist-3 amplitude (in complete analogy with twist-2 that 
(OQ)), which corresponds to longitudinally polarized virtual photon: 



A ™(Z,t) = -^J dx 



1 



H 3 (x,£,t) + ^H A {x,t,t) 
x 



k £ — x — iO £ + x — iO J |_ JV ' ' x ^ V '\ ^ 

This amplitude can be extracted from the DVCS data considering azimuthal angle de- 
pendence of the observables. 

For the new twist-3 GPDs we use WW approximation [25J 126] : 



H 3 (x,£,t) + ^H A {x,t,t) 
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where 



dH(x,£,t) ± dH(x,£,t) 



dt; dx 

After substitution of this expression into Eq. (fTHI) we obtain the expression for the twist-3 
amplitude in terms of twist-2 GPD H(x,£,t): 
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To compute the twist-3 amplitude we employ formal series decomposition of twist- 
2 GPD H(x,£,t), which represents GPDs as the infinite sum of t-channel resonance 
exchanges [28] : 
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It is easy to see that 
d±H(y,£,t) = 



dy 



(18) 



(19) 



oo n+1 



+ E E 1 



n=l,3... 2=0,2. 



x 



r z/2 



9P,(i 



Substituting this expression into Eq. (JT7I) we obtain the partial wave representation of 
the twist-3 amplitude: 

oo n+i dP(-) ^ 

n<)=nt)+2E Eg-w f (n+1)( „ +2) - po) 

n=l,3... Z=0,2... S V 1 /V 1 / 

Here A tw2 is the twist-2 DVCS amplitude given by Eqs. (l2|3l) . so that the new ingredient is 
the second term in above expression. For that term we introduce the following notation: 



oo n+1 dPl(-) 1 

A+ = 2 V V £ nJ (t)f— -ill- 77 7, (21) 

nhM... * (n+l)(n + 2)' 1 J 

and perform resummation of partial waves for it. To perform this resummation we, 
following Ref. [TJ, represent the coefficients B n i(t) as the Mellin moments of forward-like 
functions Q2v(x,t): 



B n , n +i~2v(t) = dx x n Q 2v (x,t) . 
Jo 



(22) 



Substituting this equation to ( l2Ti l changing sequence of summation and integration we 
obtain: 

MM = 2 Efe/ -^2,(M) X>' P < ( 7 > ( 23 ) 
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Q2v{%,t) with condi- 
tions 52i/ (1, t) = dS2l dx'^ = ®- ^ ^ S eas y ^° see ^ na ^ a solution of the equation has the form 



where introduced functions satisfy following equation x — 

— = 0. It is easy to see that a solutio 

of Mellin convolution: 
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Eventually, performing the summation over I in Eq. ( 1231) we obtain the following expression 
for amplitude t): 
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Here in analogy to GPD-quintessence function N(x, t) (JOj) we introduced its counterpart 
for twist-3 GPDs: 



OO .[ , OO 

n s (x, t) = j2 x2u S ^ t)= — f 1 - -) E %2v *) • ( 25 ) 

We see that at twist-3 level the function Ns(x,t) provides us with new additional infor- 
mation about forward-like functions Q2v{x,t). Therefore it is important to perform Abel 
tomography for the twist-3 DVCS amplitude. 

Abel tomography and dispersion relations for twist-3 DVCS amplitude 

Now we can obtain the function N s (x,t) if we know the twist-3 DVCS amplitude. We 
remind that the twist-3 DVCS amplitude given by Eq. ([TBI in WW approximation can 
be represented as: 



A t ^(U) = A tw2 (U) + A + (U), (26) 

where twist-2 amplitude v4 tw2 (£,t) is given by Eq. ( I2f3l) in terms of GPD-quintessence 
function N(x,t) and the amplitude A + (£,t) by the Eq. ([241 ) in terms of function Ns(x,t) 
( 1251 ). Eq. (f24j) can be brought to Abel type integral equation with help of Joukovski 
transformation and one can apply the Abel tomography described in the first section to 
the twist-3 case. The result is the following: 



N s (x,t) = -^==4 f 4 A ( WtW2 (^) - I^ tw3 (^)} • (27) 

Measuring twist-2 and twist-3 DVCS amplitudes (they can be separated studying the 
azimuthal angle dependence of various DVCS asymmetries) we can access the function 
Ng(x,t) which provides additional information about forward-like functions and GPDs. 
Using the expression for lmA tw2 ([2]) in terms of the quintessence function N(x, t) we can 
rewrite the Eq. ( [27] ) in the following form: 

m i a 1 ~ x2 ( f 1 *v 1 f 1 % ^ tw3 (^A 

It was shown in Ref. pll| that the dual parameterization representation for the twist-2 
DVCS amplitude fl2l3l) automatically satisfy the once subtracted dispersion relation: 

ReA tw2 (£,t) = 2D(t) + I vp^d( ImA tw2 (C,t) (-^ - , (29) 
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with the subtraction constant given by the D-form factor This result was obtained 
recently in Refs. [EH [H| by independent methods. We see that the dual parameterization 
automatically ensures the dispersion relations for the amplitudes. Now we can repeat the 
derivation of the dispersion relations for the twist-3 amplitude. To finish this, one has to 
substitute the expression for N s (x,t) ([271) into Eq. (1231) for the twist-3 real part of the 
DVCS amplitude. After simple calculations we obtain once subtracted dispersion relation 
for the twist-3 DVCS amplitude: 

ReA tw3 (£,t) = 2D(t) + 1 vp J d( lmA^\(,t) (-^ - ^\ . (30) 

We see that for the twist-3 amplitude the subtraction constant is given by the D-form 
factor, the same as for the twist-2 amplitude. 

Conclusions 

We extended the dual parameterization to the case of the twist-3 contributions to the 
DVCS amplitude. It is demonstrated that the twist-3 DVCS amplitude in the Wandzura- 
Wilczek approximation is determined by single function Ns(x, t)-the twist-3 analogue of 
the GPD quintessence function N(x, £)-which gives complementary information about 
forward -like functions and GPDs. We proved dispersion relations for the twist-3 am- 
plitude. It turned out that the subtraction constant for this case is the same as for the 
twist-2 case, and given by the D - form factor. 

We presented equation which directly and unambiguously express GPD quintessence 
functions N(x,t) and Ng(x,t) in terms of twist-2 and twist-3 DVCS amplitudes respec- 
tively. Relation of these equations to the Abel transform tomography is elaborated. 
Amplitudes of hard exclusive processes can be viewed as "photographs" of some object 

with axial symmetry in plane of variables X = J — \ and Y = It remains to 

decipher the deep meaning of these variables and corresponding axial symmetry. 
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